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ABSTRACT 

Aims. The goal is to determine the composition of Pluto's atmosphere and to constrain the nature of surface-atmosphere interactions. 
Methods. We perfoim high-resolution spectroscopic observations in the 2.33-2.36 /jm range, using CRIRES at the VLT. 
Results. We obtain (i) the first detection of gaseous methane in this spectral range, through lines of the + V4 and vi + V4 bands (ii) 
strong evidence (6-cr confidence) for gaseous CO in Pluto. For an isothermal atmosphere at 90 K, the CH4 and CO column densities are 
0.75 and 0.07 cm-am, within factors of 2 and 3, respectively. Using a physically-based theimal structure model of Pluto's atmosphere 
also satisfying constraints from stellar occultations, we infer CH4 and CO mixing ratios qcHj = 0.6^^ ^ % (consistent with results from 
the 1.66 /jm range) and qco = 0.5^^25^10 The CO atmospheric abundance is consistent with its surface abundance. As for Triton, 
it is probably controlled by a thin, CO-rich, detailed balancing layer resulting from seasonal transport and/or atmospheric escape. 
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1. Introduction 

Along with Triton's, Pluto's tenuous atmosphere is a benchmark 
example of an atmosphere controlled by vapor pressure equilib- 
rium with surface ices. Pluto's surface is dominated by N2 ice, 
but other detected surface compounds include CH4, CO and, ten- 
tatively, C2H6 and other non-volatile species (e.g. Cruikshank et 
al. 1997, Doute et al. 1999, Grundy and Buie 2001, 2002, Olkin 
et al. 2007, Protopapa et al. 2008, Meriin et al. 2010). Pluto's at- 
mospheric composition remains poorly known. Based on its sur- 
face abundance and large volatility, and albeit not observed spec- 
troscopically, N2 is known to dominate Pluto's atmosphere. The 
only other gas phase species to have been detected is methane. 
It was first observed in its 2v3 band near 1.66 yum by Young et 
al. (1997), providing the first estimate of its CH4 column density 
(1 .2 cm-am within a factor of 3-4, assuming a gas temperature of 
100 K). Much improved observations obtained with the CRIRES 
infrared echelle spectrograph at the VLT in 2008 (Lellouch et al. 
2009), including the detection of numerous lines, up to J = 8, 
made it possible to separate temperature and abundance effects 
in the Pluto spectra, and to constrain the depth of Pluto's atmo- 
sphere. Assuming an isothermal atmosphere, the data indicated 
a mean gas temperature T = 90^jg K and a methane column den- 
sity acHi = 0-'^5;^Q3Q cm-am. Although Pluto's surface pressure is 
uncertain, the combined analysis of these data with stellar occul- 
tation curves indicated a CH4/N2 mixing ratio qc/^j =0.5+0.1 %. 

The other species to be expected in Pluto's atmosphere 
based on volatility considerations is carbon monoxide. Searches 
for CO have been conducted akeady a decade ago (Bockelee- 
Morvan et al., 2001, Young et al. 2001), but as detailed below, 
results have been relatively unconstraining. More recently (July 
2009), and using again CRIRESA'LT, Lellouch et al. (2010) 
detected CO in Triton's atmosphere from its CO(2-0) band at 



2.35 jum, with a column density of ~0.30 cm-am, coiTespond- 
ing to CO/N2 ~ 6x10"'* for an estimated 40 /ibar surface pres- 
sure. They also searched for CO on Pluto in the (3-0) band near 
1.57 fim. The choice of that band was motivated by the much 
larger brightness of Pluto in H vs K band (geometric albedo 
~0.6 at 1.57 jum vs -0.2 at 2.35 fim). On the other hand, lines 
of the CO(3-0) band, whose strengths are ~2 orders of magni- 
tude weaker than those of the (2-0) band, are expected to show 
up only for large CO abundances, and as a matter of fact, only 
an upper limit of 1 cm-am was obtained, i.e. CO/N2 < 5x10""^ 
for a typical 15 /ibar pressure. Based on considerations on the 
surface/atmosphere equilibrium and the CO surface abundance, 
estimates of the CO abundance on Pluto are ~ 10 times smaller. 
Reaching this sensitivity in the near-IR demands the use of 
CO(2-0) band, but given that Pluto is ~3.5 times fainter than 
Triton at 2.35 fim, this in turn requires a deep integration. We 
here report on such observations. 



2. VLT/CRIRES observations 

Spectroscopic observations of Pluto were obtained during two 
half-nights on July 27 (0h40-4h25 UT) and 29 (0hl0-5hl6 
UT) 2010, using the CRIRES instiTiment (Kaufl et al. 2004) 
installed on ESO VLT (European Southern Observatory Very 
Large Telescope) UTl (Antu) 8.2 m telescope. We used the in- 
strument in adaptive optics mode and with a slit of 0.4". The 
spectral resolution, measured from the linewidths and the source 
effective size, was 60,000. We targeted a portion of the (2-0) 
band of CO, also encompassing a fraction of the vi + v\ and 
V3 + V4 bands of CH4. The instrument includes four detectors. 
On July 27, we covered the 2312-2325, 2330-2342, 2345-2356 
and 2359-2370 nm ranges. However, detectors 1 and 4 were af- 
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Fig. 1. Black: Pluto spectrum at 2330-2342 and 2345-2356 run (arbitrary units). The total integration time is 7h20 min. The spectral 
resolution is 60,000. Red, blue and green curves show synthetic spectra for an isothermal Pluto atmosphere at T = 90 K, and 
including telluric and solar lines. Green: no methane and CO in Pluto's atmosphere. Blue: methane is included with a 0.75 cm-am 
column density. Red: CO is also included with a column density of 0.07 cm-am. Blue dots locate the methane features that were 
used in the line coaddition process shown in Fig. 2. 



fected by flux losses along the slit. On July 29, we used a fast 
("windowed") readout mode in which only detectors 2 and 3 are 
read (and windowed), resulting in a noticeable gain in sensitivity 
for faint targets. Pluto's mean (East) longitude for the two dates 
was 118° and 5°, respectively. A large Doppler shift (~15.9 and 
~16.7 km/s for July 27 and 29) ensured proper separation of the 
target Pluto lines from the telluric absorptions, particularly from 
CH4. Each night we observed a telluric standard (HR 5917) to 
check for wavelength calibration and atmospheric transmission. 
All data were reduced using the standard steps of the CRJRES 
pipeline, including corrections for darks, flatfield, image recom- 
bination, replacement of bad pixels and outliers, and spectral ex- 



traction. Spectra from the two dates were finally coadded, pro- 
viding a mean S/N of 15-20 per resolving element. 



3. Modelling and evidence for CO 

The total spectrum in the 2330-2342 and 2345-2356 nm ranges 
(Fig. 1) shows the detection of many features due to the V3 + 
V4 band of CH4 (plus some lines of the vi + V4 band longwards 
of 2350 nm) representing the first detection of these bands in 
Pluto's atmosphere. As in our previous studies (Lellouch et al. 
2009, 2010) we constructed a direct, disk-integrated, line-by-line 
atmospheric model of Pluto, including solar lines reflected off' 
Pluto's surface as well as the telluric transmission and account- 
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Fig. 2. Black: Coadded CO (left) and CH4 (right) line, obtained from coaddition of the observed spectrum within +0.5 nm of (i) the 
R2-R7 and P1-P3 lines of CO and (ii) the 12 CH4 features identified by dots in Fig. 1. Colored lines: result of the same treatment 
applied to synthetic spectra. In the right panel, the feature at -0.13 nm is mostly due to the telluric counterpart of the coadded 
methane line. 



ing for the proper Doppler shifts. In a first step, we assumed 
a single, homogeneous layer at 90^^^ K, as inferred from the 
1.66 yum spectrum (Lellouch et al. 2009). This yielded a best fit 
methane column density ac/Zj = 0.75^Q g cm-am, consistent with 
the 2v3 inferences. Here the large eiTor bar combines quadrati- 
cally a factor-of-2 uncertainty due to the S/N limitations and a 
similar uncertainty due to the range in mean gas temperature. 
The lower precision compared to results from the 1.66 fim band 
is due to the lower S/N of the 2.35 fim spectrum. 

Pluto's thermal stiTicture is characterized by a "waim" 
(~100 K) upper atmosphere, a sharp inversion layer in the 3- 
10 //bar pressure range, and possibly but not necessarily a cold 
(~ 36 - 45 K) troposphere. By combining their 1.66 fim CH4 
spectrum with occultation data, Lellouch et al. (2009) found that 
the stratospheric gradient is in the range 3-15 K/km, the (cur- 
rent) surface pressure is 6.5-24 /vbar, and the troposphere is at 
most 17 km deep. Each of the allowed thermal profiles is asso- 
ciated with a methane column density acHt, required to fit the 
1.66 fim spectrum. Deeper (i.e. colder) models require larger 
methane columns than shallower models. The range for acH4 
in these models was then 0.65-1.3 cm-am, but the CH4 mixing 
ratio could be accurately determined to qc/Zj =0.5+0.1 %. We 
re-modelled the 2.35 fim spectrum with one "nominal" thermal 
profile, having a stratospheric gradient of 6 K/km, a wet tropo- 
spheric adiabat, and a surface pressure of 15 //bar, correspond- 
ing to a N2 column density of 210 cm-am (this profile is one 
of the "green" profiles shown in Fig. 4 of Lellouch et al. 2009). 
Assuming that methane is vertically uniform, an identical fit to 
that shown in Fig. 1 is then obtained for acH4 = 1.2 cm-am, i.e. 
qcHt - 0.57 %. Refitting the data with atmospheres restricted to 
p = 6.5 //bar or extending to p - 24 //bar (Table 1) indicates 
that S/N limitations dwarf thermal profile uncertainties, and the 
methane mixing ratio indicated by our spectrum is 0.6;^y 3 %. 

Besides the CH4 features, Pluto's atmosphere spectrum 
shows absorptions at the positions of most of the CO(2-0) lines 
covered by the observed spectral range (R8 — R2, RO and Pl- 
P3), and that are not fit by the previous methane-only model. 
Including CO at 90 K and with a column density of 0.07 cm- 
am clearly improves the fit in the vicinity of all these lines. 
Admittedly, several mismatches remain, that are comparable to 
or sometimes even larger than the features tentatively identi- 
fied as due to CO; examples occur near 2330.6, 2336.6, 2337 
or 2353.9 nm. To assess the significance of the result, we coad- 



Table 1. Best fit CH4 and CO abundances vs surface pressure. 



Pressure CH4 CO 

(//bar) Col. dens. Mixing ratio Col. dens. Mixing ratio 

6.5 0.6 cm-am 6.0x10"^^ 0.07 cm-am 7.0x10"'* 

15 1.2 cm-am 5.7x10"^^ 0.10 cm-am 4.8x10"'* 

24 1.7 cm-am 5.3x10"' 0.14 cm-am 4.4x10"'* 



ded the data over +0.5 nm broad intervals centered on all the 
CO line positions, without any a priori exclusion (except for the 
RO and R8 lines that are too close to detector edge). This treat- 
ment (Fig. 2, left) produces clear evidence for an absorption line 
at the "zero wavelength", around which significant structure is 
observed. Performing the same treatment on synthetic spectra 
allows us to interpret this residual structure as being due to the 
more or less random co-addition of other lines (telluric, solar, 
or Pluto CH4) in the vicinity of the CO lines. The remaining 
standard deviation between model and obser-vation (~ 1.3x10"^ 
in the units of Fig. 2, i.e. 1.8 % of the continuum) indicates a 6- 
cr detection of the "combined CO line". For an isothermal 90 K 
atmosphere, the best fit is achieved for a CO column density 
aco = 0.07 cm-am. For the above nominal thermal profile, it is 
obtained for aco = 0.10 cm-am, i.e. qco = 4.8x10""^, assuming 
uniform mixing. Given the sensitivity of the models to the CO 
abundance - and the small impact of thermal profile uncertain- 
ties (Table 1) - the allowed range is qco = 5^2*5 -^10"'*. Although 
the S/N of the Pluto spectrum is considerably lower than that 
for Triton (Lellouch et al., 2010), the accuracy on the CO abun- 
dance is comparable. This stems from the fact that the Pluto CO 
lines are instrinsically less narTow and saturated than their Triton 
counterparts, given the higher gas temperature ( ~90 K vs ~50 K) 
and ~4 times smaller CO column density on Pluto vs Triton. 

A similar line coaddition process can be performed on the 
CH4 lines. However, the larger complexity of the methane spec- 
trum implies some a priori selection of the CH4 lines. An exam- 
ple is shown in the right panel of Fig. 2, using the twelve CH4 
features marked in Fig. 1. The "combined CH4 line" is then de- 
tected at ~ll-cr. With CH4/N2 = 0.6 %, the CH4/CO ratio in 
Pluto's atmosphere is nominally equal to 12, but could be in the 
range 2^8 for the CH4 determined in this study, or 2.5-24 if the 
more precise value qcH4 = 0.5+0.1 % is used. 
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4. Discussion 

The CO amount we infer is much smaller than values previously 
reported. Using the same CO(2-0) band and IRTF/CSHELL, 
Young et al. (2001) obtained a thermal-profile dependent up- 
per limit aco < (1.2-3. 5)xlO^' cm"^, i.e. <45-130 cm-am. 
This htnit was improved by ~2 orders of magnitude from 
our VLT/CRIRES observations of the CO(3-0) band (aco < 1 
cm-am, Lellouch et al. 2010), but the present CO determina- 
tion is another factor-of-10 lower. At millimeter wavelengths. 
Bockelee-Morvan et al. (2001) obtained a tentative (formally 
4.5 cr) detection of the CO(2-l) line at 1.3 mm, using the IRAM 
30-m telescope, that, if real, would indicate a CO/N2 mixing 
ratio in the 1.2-7 % range (also thermal-profile dependent). 
However, it was more cautiously interpreted as an upper limit, 
due to strong galactic contamination in their Pluto spectrum. The 
associated CO column was <3-7 cm-am, over a factor of 6 more 
constraining than the Young et al. (2001) results. The present 
evidence for CO suggests that the authors were indeed wise to 
regard their result as an upper limit. 

Pluto's surface pressure reflects subhmation equilibrium for 
N2 ice, and N2 may follow saturation in a putative Pluto tro- 
posphere. Yet, the large CH4 mixing ratio at ~90 K in Pluto's 
stratosphere implies that it is not severely affected by atmo- 
spheric condensation, and that if Pluto has a troposphere, strong 
supersaturation of CH4 occurs (Lellouch et al., 2009). Rannou 
and Durry (2009) reached the same conclusion from microphys- 
ical arguments, but also mentioned that N2 and CO can easily 
condense. Given its equilibrium vapor pressure, only a factor of 
~8 lower than N2 at 37 K, a -0.05 % CO/N2 mixing ratio implies 
instead that CO condensation does not occur. 

As first pointed out by Lellouch (1994), CO is a cooling 
agent in Pluto's atmosphere through radiation in its pure ro- 
tational lines, but its importance has been initially overstated. 
Lellouch found that when CO cooling was included in the 
Yelle and Lunine thermal model (1989) - that included ab- 
sorption in the 3.3 jum band of CH4 as the only energy source 
- the ~100 K stratospheric temperature required a very large 
(> 10 %) CH4 abundance. Strobel et al. (1996) showed that 
this requirement was relaxed to ~ 3 % when heating in the 1 .6 
and 2.3 /im CH4 bands was considered. Still, in the presence 
of CO, a slightly negative temperature gradient (mesosphere) is 
expected in the sub-microbar region, evidence for which was 
claimed in occultation-derived thermal profiles (Elliot et al. 
2007). However, the recent analysis by Zalucha et al. (2011) 
did not confirm the need for CO to explain occultation data. 
These authors directly modelled occultation light-curves with 
physicaUy-based thermal profiles calculated from a radiative- 
conductive model inherited from Strobel et al. (1996), and in 
which the surface pressure and the CH4 and CO mixing ratios 
are free parameters. Zalucha et al. (2011) deduced CH4 mixing 
ratios of 0.18-0.94 % for occultations recorded in 1988, 2002, 
2006, and 2008 (albeit with no specific trend with time), broadly 
consistent with the observed ~0.5 % mixing ratio. In contrast, 
and although they confirmed that CO does affect the thermal 
structure in the jubar region (see e.g. their Fig. 8), their simula- 
tions indicated that occultation lightcurves are essentially insen- 
sitive to the CO mixing ratio. As a further testimony of the less 
than originally thought importance of CO for controlling Pluto's 
atmosphere state, Strobel (2008) found that CO plays at most a 
minor role in the atmospheric escape rates. 

Our inferred atmospheric CO mixing ratio is consistent with 
the long-standing predictions by Owen et al. (1993), Lellouch 
(1994) and Strobel et al. (1996). These estimates, which spanned 



the range (2-20)xl0 were obtained assuming an ideal N2- 
CO-CH4 solid solution (Raoult's law) and using the then avail- 
able CO/N2 ice mixing ratio (0.5 %, from Owen et al.). However 
a refined analysis of the near-lR spectra indicates that CO is 
present on Pluto's surface with a mixing ratio of 0.08-0.2 % 
relative to N2 (Doute et al. 1999). Using modem vapor pressure 
data (Fray et al. 2010), the ideal mixture case leads to an atmo- 
spheric CO/N2 of (1.0-2. 7)xl0""*, only marginally in agreement 
with the observed (2.5-15)xl0"'* range. Rather, the atmospheric 
and surface CO/N2 mixing ratios are the same within error bars, 
a situation predicted by the "detailed balanced" model (Trafton, 
1990; Trafton et al. 1998). In this scenario, surface-atmosphere 
exchanges in presence of escape and seasonal transport lead to 
an atmospheric composition reflecting that of the accessible ice 
reservoir from which it is replenished. In the case, relevant for 
a CO / N2 mixture, where no fractionation occurs during escape 
or transport, the process equalizes the mixing ratios in the atmo- 
sphere and the volatile reservoir. This is made possible by the 
formation of a thin CO-enriched surface veneer in equilibrium 
with the atmosphere according to Raoult's law. Lellouch et al. 
(2010) found this scenario to be valid for CO on Triton. 

The altemative scenario would be that the enhanced CO 
abundance compared to the ideal solution case results from the 
presence of pure CO patches. The CO partial pressure is 4-25 
nbar, i.e. 80-500 times less than the vapor pressure of pure CO. 
Hence, patches of pure CO covering 0.2-1.2 % of Pluto's surface 
could in theory produce the elevated CO abundance. Problems 
with this scenario are that (i) the formation of pure CO grains 
is not expected thermodynamically, given the complete misci- 
bility of N2 and CO (ii) CO is not buoyant in N2, inhibiting the 
sublimation of any CO patch (Stansberry et al. 1996). Yet, the 
longitudinal distribution of the CO 1.58 fim ice features shows 
evidence for a CO-rich region near longitude L = 1 80° (Grundy 
and Buie 2001), which (coincidentally or not) corresponds to the 
brightest and least red region in HST maps near Pluto's equator. 
(Buie et al. 2010). This situation contrasts with the Triton case, 
where the longitudinal distributions of the N2 and CO ice bands 
are remarkably similar (Grundy et al. 2010). How could CO be- 
come concentrated in a locaUzed spot on Pluto remains to be 
elucidated. Hopefully, combined imaging and spectroscopy from 
New Horizons will reveal the nature of this CO-rich region and 
whether it may contribute to the CO atmospheric abundance. 
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